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ABSTRACT 

Likely progenitors for the GRBs (gamma-ray bursts) are the mergers of 
compact objects or the explosions of massive stars. These two cases have 
distinctive environments for the GRB afterglow: the compact object explosions 
occur in the ISM (interstellar medium) and those of massive stars occur in the 
preburst stellar wind. We calculate the expected afterglow for a burst in a 
Wolf-Rayet star wind and compare the results to those for constant, interstellar 
density. The optical afterglow for the wind case is generally expected to decline 
more steeply than in the constant density case, but this effect may be masked 
by variations in electron spectral index, and the two cases have the same 
evolution in the cooling regime. Observations of the concurrent radio and 
optical/X-ray evolution are especially useful for distinguishing between the two 
cases. The different rates of decline of the optical and X-ray afterglows of GRB 
990123 suggest constant density interaction for this case. We have previously 
found strong evidence for wind interaction in SN 1998bw/GRB 980425 and 
here present a wind model for GRB 980519. We thus suggest that there are 
both wind type GRB afterglows with massive star progenitors and ISM type 
afterglows with compact binary star progenitors. The wind type bursts are 
likely to be accompanied by a supernova, but not the ISM type. 



Subject headings: gamma-rays: bursts — stars: mass loss — stars: supernovae: 
general 



1. INTRODUCTION 



Although the study of GRBs (gamma-ray bursts) has been revolutionized in the 
past few years by finding a number of precise positions and distances, the nature of their 
progenitor objects remains uncertain (see Meszaros 1999 for a review). The production of a 



- 2 - 



large amount of energy in a short time has naturally led to models involving compact objects 
(neutron stars and black holes). Fryer, Woosley, & Hartmann (1999) have summarized 
possible progenitors involving black hole accretion disks: neutron star - neutron star binary 
mergers (NS/NS), black hole - neutron star mergers (BH/NS), black hole - white dwarf 
mergers (BH/WD), massive star core collapses, and black hole - helium star mergers 
(BH/Hc). Based on estimated formation rates and on accretion disk models with high 
viscous forces, Fryer et al. (1999) suggest that NS/NS and BH/NS mergers dominate the 
population of short- duration GRBs and that massive stars and BH/He mergers dominate 
the long-duration bursts. The afterglows observed to date would have massive star 
progenitors in this scenario because they followed long-duration bursts. 

One way of distinguishing between the progenitor models is to examine the position of 
the GRB in the parent galaxy. Paczyhski (1998) pointed out that NS/NS binaries would 
be expected to have a significant space velocity, which would carry them many kpc from 
their birthplaces. The observational evidence for the association of several GRBs with 
star forming regions then provided weak evidence against the NS/NS merger progenitors 
and favored massive star progenitors. The population synthesis calculations of Fryer et al. 
(1999) supported this conclusion. However, Bloom, Sigurdsson, & Pols (1999a) found a 
time to NS/NS merger of ~ 10 8 years. These objects would then follow the star formation 
rate, although ~ 15% of them might occur well outside of dwarf galaxy hosts. 

The connection of GRBs to massive stars became stronger with the discovery of the 
Type Ic supernova SN 1998bw in the error box of GRB 980425 (Galama et al. 1998). The 
high energy inferred for the optical supernova, (2 — 3) x 10 52 ergs (Iwamoto et al. 1998; 
Woosley, Eastman, & Schmidt 1999), and the high expansion velocity inferred for the 
radio supernova (Kulkarni et al. 1998) strengthen the GRB connection. Li & Chevalier 
(1999) found that the evolution of the radio source indicated non-uniform energy input to 
the blast wave, as is also inferred in GRBs. They also found evidence that the radio SN 
1998bw interacted with the stellar wind expected from the massive star progenitor. In this 
paper, we emphasize that a stellar wind environment is an unavoidable consequence of a 
massive star progenitor and that the nature of the GRB afterglow emission can provide 
a discriminant between massive star and compact binary progenitor models. In § 2, we 
discuss the expected afterglow for a massive star progenitor model and in § 3 place our 
models in the context of observations. Our discussion concentrates on the cases s = 2 
(stellar wind) and s = (interstellar medium), where the ambient medium has p oc r~ s . 
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2. AFTERGLOWS IN A MASSIVE STAR WIND 

In the existing massive star GRB progenitor models, the most likely progenitor 
is the stripped core of a massive star, i.e., a Wolf-Rayet star. MacFadyen & Woosley 
(1999) consider single massive stars whose cores directly collapse to black holes. The 
stars have an initial mass >25 M , which is the type of star that is likely to lose its H 
envelope in winds. Paczyhski (1998) noted that the requirement of a rapidly rotating core 
might necessitate a close binary companion, which again points to Wolf-Rayet stars. The 
winds from these stars in our Galaxy have velocities of 1,000 — 2,500 km s" 1 and mass 
loss rates M tv 10~ 5 — 10~ 4 M yr _1 (Willis 1991). On evolutionary grounds, Langer 
(1989) advocated M~6x 10~ 8 (M WR / M ) 2 - 5 M yr" 1 , where M WR is the mass of the 
Wolf-Rayet star. If the stellar mass drops to ~ 3 M at the end of its life because of mass 
loss, M ~ 10~ 6 M Q yr^ 1 at that time. The host galaxies of GRBs are likely to be of low 
metallicity. Willis (1991) notes that there is no evidence for a metallicity dependence of 
mass loss from a particular type of Wolf-Rayet star, but that metallicity may affect the 
distribution of Wolf-Rayet types. 

The effects of Wolf-Rayet winds can be observed in the wind bubbles around some 
of these objects. In some cases where the bubble is well-studied, the bubble expansion 
and X-ray emission are consistent with a weaker wind than that inferred from direct 
observations of the Wolf-Rayet star (e.g., NGC 6888, Garcia-Segura & Mac Low 1995). The 
position of the wind termination shock, r t , can be estimated by equating the ram pressure 
in the wind with the pressure in the bubble. We estimate r t /Rf, m (2Rf,/v w t w y/ 2 , where 
Rt, is the radius of the wind bubble and t w is its age. For Rf, = 3 pc, t w = 2 x 10 4 yr, and 
v w = 1, 000 km s -1 , we have r t ~ 5 x 10 18 cm. At r t , the density increases by a factor of 4 
and becomes approximately constant at larger radius. The radius r t is expected to increase 
as the bubble evolves and its pressure drops. At the end of its life, a Wolf-Rayet star is 
thus expected to be surrounded by a substantial p oc r~ 2 medium. The radio observations 
of Type Ib/c supernovae can be interpreted as interactions with such a medium in models 
with constant efficiencies of production of relativistic electrons and magnetic fields (Weiler 
et al. 1999; Chevalier 1998). 

For the wind density, we take p = Ar~ 2 , where A = M/4irv w = 5 x 10 n A* g cm -1 . The 
reference value of A corresponds to M — 1 x 10~ 5 M yr -1 and v w = 1000 km s -1 . If a burst 
were able to occur in a red supergiant star, the low wind velocity, v w ~ 10 km s -1 , would 
lead to a higher circumstellar density because M>10~ 7 M yr~ x up to ~ 10 -4 M yr -1 
is expected. The only supernova case where we have inferred a lower circumstellar density 
is around SN 1987A, which exploded as a B3 I star and was thus too cool to drive a 
strong wind with its ultraviolet radiation field. In this case, the r~ 2 wind terminated 
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at ~ (3 — 4) x 10 17 cm and the supernova radio flux rose sharply when the shock front 
encountered denser gas (Ball et al. 1995; Chevalier 1998). 

The scaling laws that are appropriate for GRB interaction with an s = 2 medium have 
been described by Meszaros, Rees, & Wijers (1998) and by Panaitescu, Meszaros, & Rees 
(1998). Our aim here is to examine the specific predictions for interaction with a Wolf-Rayet 
star wind. We calculate the expected emission for a thin shell model (cf. Li & Chevalier 
1999). We adopt a nucleon-to-electron number density ratio of two, which is appropriate 
for winds of Wolf-Rayet stars that are predominantly helium (and perhaps carbon/oxygen). 
For the mass loss case, a density of 1.67 x 10~ 24 g cm~ 3 is attained at r pa 5.5 x 10 17 Ay 2 
cm. For an adiabatic blast wave in an s = 2 medium (Blandford & McKee 1976), we find 
7 2 = R/Act and R = {9Et/A7iAc)^ 2 = 2.0 x 10 17 E^A^ 2 ^ cm, where 7 is the Lorentz 
factor of the gas, R is the observed radius near the line of sight, E 52 = E/10 52 ergs is the 
explosion energy, and t = t day day is the time in the observer's frame. Over the typical time 
of observation of a GRB afterglow, the shock front is expected to be within the s = 2 wind, 
except for the unusual case of a progenitor like that of SN 1987A. 

In our model, the synchrotron emission frequency of the lowest energy electrons is 

v m = 5 x 10 12 (^) V2 (^e/0.1) W0.1) 1/2 i4%af Hz, (1) 
and the flux at this frequency is 

F - = 20 ( ^iT (^ £ ) 1/2 ( £ "/°- 1 ) ,/2B - 2 ^C /2 »uy. (2) 

where z is the redshift, and e e and are the electron and magnetic postshock energy 
fractions. These expressions assume a flat universe with H Q = 65 km s" 1 Mpc -1 . As in 
discussions of afterglows in the ISM, the magnetic field is assumed to be amplified by 
processes in the shocked region and is not directly determined by the field in the Wolf-Rayet 
star wind. For electrons with a power law energy distribution above r y em , dN e /d^ e oc r y~ p , 
the flux above v m is F v = F Vm (y jv m )-^l 2 oc e P-i e g+ 1 )/ 4 f;(p+i)/4^ r (3p-i)/4 I/ -(p-i)/2 j 

provided the electrons are not radiative. F u = F Um {v / v m ) 1 ^ below v m until the spectrum 
turns over due to synchrotron self-absorption at 

va - 1 x 10 11 (^)" /5 {^AY\e B I^E^ATt2!; Hz. (3) 

We estimate the effects of synchrotron cooling following the discussions of Sari, Piran, & 
Narayan (1998) and Wijers & Galama (1999) for the s = case. Radiative cooling becomes 
important at a frequency 

u c « 1 x 10 12 (^)~ 3/2 (es/OAr^ElfA: 2 ^ Hz, (4) 
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provided v c > u m . For v > v c , F v oc e v-\%-^^E^)lH-^~^l A v^l 2 . 

These results can be compared to similar results for interaction with a constant density 
interstellar medium with n Q m I cur 3 (e.g., Waxman 1997; Sari et al. 1998). In that case 
F Um is constant with time, v m is comparable to the value given above, va is independent of 
time, and v c decreases with time. If va, v m , and v c have the same relation to each other in 
the s = and s = 2 cases, the appearance of the spectrum at one time is similar for both 
cases, but the evolution is different. At high frequency (optical and X-ray), for s = the 
flux evolution goes from adiabatic (oc £~( 3 p~ 3 )/ 4 ) to cooling (oc £~( 3 p~ 2 )/ 4 ) while for s = 2 the 
flux evolution goes from cooling (oc £~( 3 p~ 2 )/ 4 ) to adiabatic (oc £~( 3 p _1 )/ 4 ). While cooling, the 
two cases have the same spectrum and flux decline. At low frequency (radio) with v < u m , 
the flux evolves as t 1 / 2 for s = 0, but can make a transition from oc t to constant for s = 2. 

3. COMPARISON WITH OBSERVATIONS 

We have noted the evidence that SN 1998bw is interacting with a circumstellar 
wind (Li & Chevalier 1999), and here discuss the cosmological GRBs. The more detailed 
models that have been developed for comparison with the best observed GRB afterglows 
have taken a low, constant density surrounding medium as the starting point. For GRB 
970508, Wijers & Galama (1999) and Granot, Piran, & Sari (1999) found n Q = 0.030,5.3 
cm- 3 , E 52 = 3.5 , 0.53, e e = 0.12,0.57, and e B = 0.089,0.0082, respectively. The difference 
between these results reflects, in part, the difficulty of accurately determining the significant 
frequencies, such as va- The low density found for this case is indicative of an interstellar 
density. However, these models primarily depend on the spectrum at one time, which can 
be the same for the s = and s = 2 cases. Panaitescu et al. (1998) have presented a 
successful s = model for the evolution of GRB 970508, but a number of deviations from 
the simple model are included, so the case is not clear. 

With s — 0, the optically thin, adiabatic flux evolution can be described by F v oc t a v^ 
with a = 1.5(3. This type of evolution was observed in the afterglow of GRB 970228 
(Wijers, Rees, & Meszaros 1997). For s — 2, the expected power law evolution has 
a = (3/9 — l)/2 for a constant energy blast wave, so a mechanism must be found to flatten 
the time evolution if this GRB evolved in a circumstellar wind. The effects of beaming and 
a change to nonrelativistic flow steepen the time evolution, but continued power input from 
ejecta can flatten it. This effect is inferred to occur in radio supernovae (Chevalier 1998). 
Following Rees & Meszaros (1998) for the s = case with power input from ejecta with a 
mass gradient M(>Tf) oc Tj n where Tj is the Lorentz factor in the freely expanding ejecta, 
the evolution follows F v oc t~[ 2_ ' 3 ( n + 6 )]/( Tt + 4 ) for the s = 2 case. For (3 = —0.75, the property 
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a = 1.5/9 is recovered when n = 5.33. The optically thin, v > v m evolution would then be 
the same as in the constant density (s = 0) case, but other aspects of the evolution would 
be different. In particular, the expansion would be much less decelerated and the apparent 
radius of the blast r oc t {n+2)/{n+A) oc t a79 and its Lorentz factor 7 oc r 1/(n+4) oc r an . At 
optically thick wavelengths, the flux would increase as F u oc t L57 as opposed to F u oc t 1 ^ 2 
for the constant density case. Radio observations could distinguish this kind of evolution, 
but GRB 970228 was not detected in the radio. 

Kulkarni et al. (1999) noted that the afterglow of GRB 990123 was consistent with 
p = 2.44 in a s = model, where p is the electron energy spectral index. The optical 
emission declined with a = —1.10 ± 0.03, suggesting adiabatic evolution, while the X-rays 
declined with a = —1.44 ± 0.07, suggesting cooling evolution. The faster decline at X-ray 
wavelengths is distinctive of s = evolution and is the opposite of expectations for s = 2. 
Most of the well observed optical afterglows also have a in the range -1.1 to -1.3, which 
is plausibly modeled by blast wave evolution in a constant density medium. However, 
the possibilities of flat electron spectra or cooling evolution do not allow s = 2 models to 
generally be ruled out. 

An afterglow with a steep decline was that of GRB 980519. Its optical emission 
followed F u oc t-( 2 - 05±0 - 04 ) ; which, with the observed (3 = —1.05 ± 0.10 for optical and X-ray 
data, is consistent with expansion in an s = 2 medium (Halpern et al. 1999). The steep 
decline might be the result of beaming (Halpern et al. 1999; Sari, Piran, & Halpern 1999) 
instead of interaction with a wind. One way of distinguishing the wind interaction case is 
to observe the optically thick radio flux rise, F u oc t for a wind as opposed to F v oc t 1 / 2 for 
a constant density. The effect of beaming in the s = case is to decrease the expansion 
and thus to increase the difference between the models. Radio observations of GRB 980519 
were made at ages 0.3, 1.1, and 2.8 days at 8.3 GHz (Frail et al. 1998). Fig. 1 shows that 
an s = 2 model is capable of fitting the radio, as well as the optical and X-ray, data. 

The model light curves in Fig. 1 are obtained using the synchrotron self-absorption 
model of Li & Chevalier (1999). The model takes into account the dynamical evolution of a 
spherical, constant-energy blast wave in a pre-burst s = 2 stellar wind, relativistic effects 
on radiation, and synchrotron self-absorption. The distance to GRB 980519 is unknown, 
so we take a fiducial value of z = 1. Adopting a relatively low magnetic energy fraction 
of eb = 10~ 5 , we find that the following combination of parameters fits all available data 
reasonably well: E 52 = 0.54, e e = 0.62, p = 3, and A* = 4.3. Rough scalings of these 
parameters for other choices of €3 are given in Li & Chevalier (1999). The value of p is 
higher than that normally found in GRB afterglows, but is within the range found in radio 
supernovae (Chevalier 1998 and references therein). The model shown in Fig. 1 does not 
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include synchrotron cooling. The good fit to the X-ray data is then expected because the 
spectral index joining optical and X-ray emission, (5 = —1.05 ± 0.10 (Halpern et al. 1999), 
is consistent with the separate optical and X-ray indices. From eq. (4) and the parameters 
for GRB 980519 given above, we have v c 4 x 10 16 Hz for tday = 1, which is below the 
X-ray frequency (3 keV = 7 x 10 17 Hz). We consider the agreement adequate in view of the 
expected gradual turnover in flux and the observational and theoretical uncertainty in the 
value of v c . The model does require a low value of e B) as has also been inferred in a number 
of other afterglows (Galama et al. 1999). 

The model R-band flux densities on days 60 and 66 of Fig. 1 fall short of the observed 
values by nearly two orders of magnitude; the emission on these two days is presumed to 
come from the host galaxy of the burst (Sokolov et al. 1998; Bloom et al. 1998a). However, 
the observed flux densities on days 60 and 66 are close to those expected of SN 1998bw at a 
cosmological distance of z — 1 (see Fig. 1). Despite excellent seeing conditions at Keck II, 
Bloom et al. (1998a) found little evidence for the extension expected of a host galaxy, so 
the presence of a SN 1998bw-type supernova is possible. Observations of the source at later 
times than have been reported should be able to confirm or reject this possibility. 

Another burst with a relatively steep decline with time is GRB 980326, in which the 
optical F v oc t- 2 - 10±0 - 13 (Groot et al. 1998). There is no radio data for this object, but the 
last optical observation is a factor ~ 10 above the power law decline; this is not the host 
galaxy because it is not present at a later time (Bloom & Kulkarni 1998). A possible source 
for the emission is a supernova (Bloom et al. 1999b), as in the case of GRB 980519. Such an 
event would be consistent with expansion into a wind and the explosion of a massive star. 
Bloom et al. (1999b) found that the spectrum of the late time source is consistent with 
that of a supernova, which rules out the possibility of a rise in the nonthermal afterglow 
emission as observed in GRB 970508 (Panaitescu et al. 1998) and the radio emission from 
SN 1998bw (Kulkarni et al. 1998; Li & Chevalier 1999). 

Based on the discovery of the SN 1998bw/GRB 980425 association, Bloom et al. 
(1998b) proposed a subclass of GRBs produced by supernovae (S-GRBs). We also place 
SN 1998bw in a supernova class of GRB, but with different properties from those of Bloom 
et al. (1998b). By supernova, we mean an observed event with an optical light curve and 
spectrum similar to that of a Type I or Type II supernova. In our picture, the GRB 990123 
is an example of a GRB interacting with the ISM (interstellar medium). The similar (slow) 
rates of decline observed in other GRB afterglows suggests that this type of object is the 
most common. These objects are not interacting with winds, do not have massive star 
progenitors, and are not accompanied by supernovae. Plausible progenitors are the mergers 
of compact objects. The GRBs in the wind type are interacting with winds, have massive 
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star progenitors, and are accompanied by supernovae. SN 1998bw/GRB 980425 is the best 
example of this class of object. GRB 980519 and GRB 980326 are other possible members. 
Bloom et al. (1998b) propose that S-GRBs have no long-lived X-ray afterglow because of 
synchrotron losses. We suggest that X-ray afterglows are possible, although with steeper 
time evolution than in the ISM case, provided e# is small, as is also inferred for some of the 
ISM type afterglows. Bloom et al. further propose that S-GRBs have single pulse GRB 
profiles. We have found evidence for nonuniform energy input in SN 1998bw/GRB 980425 
(Li & Chevalier 1999) and suggest that both types of bursts are powered by central engines 
with irregular power output. The GRB itself may not allow a classification of the event; 
the two cosmological GRBs mentioned as possible members of the wind class have multiple 
peak time structure (in 't Zand et al. 1999; Groot et al. 1998). 



We are grateful to J. Bloom, S. Kulkarni, and an anonymous referee for useful 
comments and information. Support for this work was provided in part by NASA grant 
NAG5-8232. 
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Fig. 1. — Wind interaction model for the afterglow of GRB 980519. The radio, optical, 
and X-ray data are taken, respectively, from Frail et al. (1998), Halpern et al. (1999), and 
Nicastro et al. (1999). A factor-of-two increase is applied to the original optical data (open 
squares) to account for dust extinction. The X-ray flux densities (large filled squares) are the 
observed fluxes in the 2-10 kev band divided by the band width. The solid curves with small 
filled dots are the light curves from a simple synchrotron self-absorption model described in 
the text. The open circles are the R-band flux densities that SN 1998bw would have at a 
cosmological redshift of z — 1 based on its U-band flux densities at z — 0.0085 (Galama et 
al. 1998). 
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